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Abstract 

Previously developed method |l], Q for finding asymptotic solutions of Vlasov equa- 
tions using two-dimensional (in coordinate x and time t) Laplace transform is here ap- 
plied to consider transversal oscillations and waves in low-collision quasi-neutral (n, ~ 
n e ) Maxwellian electron-ion plasmas. We obtain two branches of electron waves: the 
ubiquitous one of high-frequency and high-velocity oscillations and the unusual low- 
velocity one. Taking into account Coulomb collisions in the limit m e <C rm, Vi <C v e , 
and T e m e <C Tiirii results in expressions for transversal plasma-electron oscillation/wave 
decrements with a damping of the low- velocity electron branch ~ nj^ /vt^ , where nj is 
the ion density and v e is the mean electron velocity. It ought to rehabilitate Vlasov prin- 
cipal value prescription for relevant integrals, but to supplement it with representation 
of an asymptotical solution as a sum of exponents (not a single one!). "Non-damping" 
kinematical waves in low-collision plasma transform in the damping ones at reasonably 
chosen iteration process. 

PACS numbers: 52.25 Dg; 52.35 Fp. 

Key words: plasma oscillations; plasma waves; Landau damping; electron waves; Coulomb 
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1 Introduction 



Propagation of plane transversal electromagnetic waves in plasmas is described by asymptotic 
in coordinates and time solutions of the coupled system of kinetic equations for electrons and 
ions and Maxwell equations for the electric field: 
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the plane wave is moving in x-direction, /q is Maxwell distribution function. Let us note that 
one can add to the right hand sides of Eqs. ([!]), (0) some collision integrals. 

Here we limit ourselves to a consideration of the particular case of plane electron waves 
in infinite homogenous fully ionized plasma with the boundary electric field perturbation 
E(0,t) = E exp(iut). 

In close analogy to our previous paper |I[ , where we have considered the damping of lon- 
gitudinal waves in the electron-ion low-collision plasmas, the problem is solved using iteration 
technique with substitution of collisionless distribution functions into the Coulomb collision 
integrals. 



2 The null iteration 



Neglecting the ion constituent in Eq.(^) and using Laplace transforms in Eqs.(|T])-(^) we arrive 
at 
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where af 2 = a\, 2 ± zoo and /W, F Pl , F P2 , E Pl , and E P2 are, correspondingly, Laplace 



transforms of / (e) (v, 0, t), f (e) (v, x, 0) 
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E z (Q,t), and E z (x,0). 



Neglecting for simplicity initial and boundary values 

f(i) f(i) F F E 
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(they do not affect characteristic frequencies uj = —ipi and wave numbers k = —ip 2 ) one 
obtains the following equation for the double poles in pi and p 2 . 
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where lol = j4:TTe 2 n e /m e is Langmuir frequency. Using transformation 
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where v 2 can be approximated by the mean square velocity defined by Maxwell exponent 
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one obtains 
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This implies the pole in the complex pi plane: 

Pi = uo 

and the corresponding pole in the complex p 2 plane defined from Eq.(|T5|). 

Taking account for v 2 -C c 2 one obtains from Eqs.flTB]) and (|17|) two solutions: 
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The solution (|T9|) at cj^ < cu, that is 
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is the well-known result of non-damping transversal electromagnetic high-frequency waves in 
fully ionized plasma || . Phase velocity of this mode is greater than c, but this result is by no 
means related to the applicability of Maxwell distribution function fff^ at high velocities v up 
to c, but, instead, is due only to the Maxwell field equation Eq.(f|). 



The solution (|20|) for non-damping low-velocity waves appears to be more intriguing: 
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with phase and group velocities 
which are not dependent on ojl. 
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Since we assume definiteness and convergence of the inverse Laplace transformation we 
should discuss an appearance of numerous poles (p% + p2V x ) = in integrals in dv x in Eq.(|12"D 

(1 2) 

depending on running values pi, p 2 - The results of calculations of poles p\ = iu> and p\ , 
Eqs.(|T^)-(pOD, show that at these values of p± and P2 the integral in dv x in Eqs. (|12D— (pT3|) 



appears as logarithmically divergent. Strictly speaking, approximation (|T3|) implies that this 
integral is defined in a principal-value sense. In this case the inverse Laplace transformation 
from approximate "image" to "original" function is also definite both for E piP2 and /IJL- 
The criterium for the validity of some found solutions fi(v x ,x,t) and E(x,t) is the fulfilment 
of original Eqs. ([!])- (|3|) on these functions. Such definite solutions can be obtained if one 
takes the principal value prescription for integrals in dv x (using approximation (|l"3D or not). 
This prescription is not an arbitrary agreement, but is the necessary consequence of Laplace 
transform existence. So, incorrectness of solving Vlasov equations by Vlasov himself (in case 
of longitudinal plasma oscillations [[|) was in representing the solution in the form of a single 
exponent exp(iut — ikx) rather than in the form of / dv x in principal value prescription. 
Asymptotical solution in general case is some complex function which can be expanded in a 
series of exponents, but it is not necessarily a single exponent. 
Substitution of expression 
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into Eq.(|l5|) confirms the existence of two opposite in signs poles of the "image" E piP2 . Cal- 

(2) 

culation of residua in these poles p 2 and amplitudes of the electric field oscillations for this 
low-velocity non-damping mode results here in a trivial asymptotic solution in the form of a 
standing wave 

#0M)asymp = E e tujt cos(kx) . 

However real existence and amplitudes of this mode must be defined at accounting for 
additive constituents from all other partially coupledf] or independent boundary and initial 



conditions (|TT|) which have been omitted till now for simplicity. So, one of such additional 
conditions might be, for example, even if partly absence of backward wave, etc. 

If the boundary and initial conditions are nevertheless such that some oscillatory mode 
with frequency u is represented only with a single forward wave (and/or backward wave), all 
contradictions will be removed if the logarithmically divergent integral in the partial dispersion 
relation of this mode ( (|24D or any other) will be treated in the Vlasov sense of the principal 
value. In general case including collision plasma and longitudinal waves one can get presence of 
a single travelling non-damping or damping forward wave only if the latter has non-exponential 
form F(ut — kx) with a non-exponential boundary condition, for example, E = E Q cos(ujt) 
(that is at least a sum of two complex- valued exponents). 

It might be interesting to note that quite nearly to the poles (|24[) there are located the 



values of variables P2 in the integrand E PlP2 of inverse Laplace transformation 
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with zero contribution into Laplace (^-integral in the formula for an asymptotical value of 
E(x,t). Such proximity between values of P2 variable at which E PlP2 turns to the infinity 
and to the zero is intriguing. What will happen in case when v x — > 0? One can suppose, 
for example, that this case may be accompanied with a considerable augmentation of the 
length/time needed to set up the asymptotical regime to really unobservable large values. 
Be that as it may, strict interpretation of this mode features is possible only at account for 
relativistic corrections in Eq. ([!]). 

2 One of such additional conditions might be, for example, the absence of a backward wave. 



3 Coulomb low-collision plasma 

Coulomb collision integral is included into Eqs. (|Tj) and (|1^) in the form 
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where Q P1P2 is Laplace transform of the usual Coulomb collision integral dominated by the 
electron-ion collision term 0: 
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where u = V — v and V is the ion velocity. 

Thus, the calculation of the collision contribution into E PlP2 reduces to the calculation of 
the term 
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After integration by parts and simple transformation of the type (p"3|) 
characteristic equation for determining decrement S: 
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Substituting the value p\ = ko and assuming 

P2=P? + S |*|« 

where p 2 is defined by Eq . (|i~9"D , one obtains for the coordinate damping decrement: 

S = S X - 
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with Si < for wave with k < ( PH) travelling to the right direction. 

But more interesting solution appears for the wave (|2"2]). Substituting into Eq.i 
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and keeping in Eq. (|30|) only terms to the order 5 3 , accounting for v\ <C c 2 , one obtains the 
following striking nonlinear result for the low- velocity electron branch decrement: 
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This intriguing result can be tested experimentally, but one should have in mind that in 
low-velocity/low-frequency region there might also exist other electron-ion branches defined 
by the ion current in the Maxwell equation (|||), which can complicate interpretation of test 
results. It is worth also to note that decrements Si >2 , generally speaking, should not be used 
automatically as simple additive parts of the whole collision damping decrement of partially 
ionized plasma where collisions of electrons with neutral particles are also present and should 
be taken into account. 

4 Conclusions 

Application of our method |3J of 2-dimensional Laplace transformation to plasma transversal 
oscillation equations with calculating logarithmically divergent integrals with the principal 
value prescription results in determination of plasma oscillation frequencies (or wave num- 
bers). We have obtained dispersion equations for non-damping oscillatory modes k(u) of 
quasi-neutral Maxwellian collisionless fully ionized plasma and damping modes of low-collision 
plasma including high-velocity (with phase velocity > c) and low-velocity (with phase and 
group velocities ~ w e C c) transversal modes. 

With the help of the same method we have also obtained the damping decrements of these 
modes due to Coulomb electron-ion collisions in the low-collision fully ionized plasmas. The 
most striking thing is a non-linear nature of the decrement for the low-velocity mode 



where v 2 is the mean-square velocity of the electrons. 

The obtained results on propagation and damping of plasma waves and oscillations can 
be useful not only in applications to laboratory plasma studies, but mainly in theoretical 
evaluations of non-thermal energy and damping lengths in the solar atmosphere as well as in 
interplanetary and interstellar media. 

Let us emphasize once more: there is no necessity to appeal to Landau's rule of passing 
around poles in calculations of indefinitely (logarithmically) divergent integrals both in the 
case of transverse as well as of longitudinal plasma waves. The effect of dissipative "Landau 
damping" , see Jl[ , does not exist in nature and is no more than some abstract great fiction of 
theorists. 

It ought to rehabilitate Vlasov principal value prescription of his relevant logarithmically 
divergent integrals, however to generalize his solution with using not a single exponent, 
but some combination of exponents as for an asymptotical solution as well as for functions in 
boundary and initial conditions (for example, these conditions have to be specially selected in 
order to avoid unphysical divergent at x — > oo backward waves, etc.). 

Asymptotical solution is some linear combination (a) of exponents (Laplace expansion) 
which must satisfy linear plasma equations. The linear combination (b) of the exponents 
from boundary and initial conditions ought to be considered as a selector for the exponents 
to be included in group (a). It should be noted here that the boundary condition of the type 
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E ~ Eo exp (iut) is some mathematical abstraction and can not be realized in a real physical 
situation. 

According to formulas given in this paper for collisionless plasma one has 



P2 E P1 + F P1 - - X 



Pl E P2 + F P2 _ x r (giffl ~ ^ffl) ^*M 3 

C 2 7 pi + MxP2 

with A = 47r|e|n e /c 2 and 



(37) 



2 Pi , ^lPi f d fo{e) d 3 u 



The point pi + u x pi = is not a point of singularity of integrals in Eqs.([37]) and (p8|) because 
the integrals are defined with the help of the principal value prescription. But the point 
Pi + v xP2 = is a pole of fj-^ in Eq.(^) with functions fj£, E Pl and others, which, in 
its turn, can have singularities in pi, p 2 . 

According to || §] this pole corresponds to solutions 
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of equation 
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with arbitrary (or a;), C(-y), and t^. 

Inverse Laplace transform ip(v,x,t) of expression 
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with some arbitrary initial (<£> P2 ) and boundary (y pi ) conditions gives the general complete 
solution of Eq. (f41~|) , including as well the alternative partial solutions ( |39"D and (pEOl) . So, one 
can always select from the set of solutions of Eq. ([J) 

f[ e \v, x, t; ip) = f[ e \v, x, t) + <p{v, x, t) (43) 

with arbitrary initial and boundary conditions a variant with f^> = fjg = 0. This solution, 
due to Eq . (|37D , induces appearance of electrical field E(x, t) and of unphysical backward waves. 
Let us assume the boundary electric field to be of the following form 

E(0,t) = E cQs{ut) . (44) 

Then it is natural to expect 

f[ e \v, 0, t) = a^v) cos(wt) + b u {y) sin(c^) (45) 
with some arbitrary coefficients a w (v) and b^v). 



In this case Laplace transforms are 
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According to Eq.(|37j) with account for collision damping (through the additional collision 
term in G) the term p 2 E pi /G leads to asymptotical waves with exponents of (hot ± ikx ± Sx) 
and (—icot±ikx^fSx), and the term with fj^ leads to waves with exponents of (itut±ikx±Sx) 
(with coefficient a u in Eq. (|47|) ) and (—hot ± ikx =F (with coefficient b u in Eq. (|47|) ). It is 
evident that it is always possible to eliminate backward waves by the proper selection of 
and b^. Possible more convenient modification of Eq.(O) can be also, for instance, 
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and /^p are all proportional to i?o- 
At account for the simultaneous presence of the other mode fl22|), it is seen that the already 
defined above a and b have to lead to removal of backward waves in this mode too, due to its 
analogous structure. 

The paradox consists in the fact that at damping electrical fieldft due to collision terms in 
G, there exist non-damping "ballistic" kinematical oscillations in fi (v , x, t) (in [^, Q there 
was implied, instead, collisionless Landau damping with the same result). These "ballistic" 
oscillations appear to be unremovable. But there is no mystery. If collision term, even if 
partly, included some additional term of the form 
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then, instead of denominator pi + v x p 2 and p% + u x p 2 in Eqs.( ^6|)-( |38D, one would have 
denominator p 1 + v x p 2 + C(v) or p\ + u x p 2 + C(u) and the "ballistic" term would be damping 
too in the same extent as the oscillations of electric field. 

It appears that, for instance, as the first collision iteration one could more reasonably use 
the substitution in Eq.(p6|) 
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is 



where f/ p J P2 ) is the collisionless approximation of fj^ p2 according to Eq. (|36|) , and Qp iP2 
the function of pi, p 2 , and v defined by the Laplace transformation of expression (|28|) with 
substitution f^ p2 — > (/ P j P2 ) • Here arises also the need for some new correction of coefficients 

a and b in the boundary and initial conditions fj£> ad fjg in order to avoid growing backwards 
waves. 

By this way non-damping of ballistic mode in low-collision plasma can be related with some 
approximation made in iteration procedures (yet nobody has still investigated convergence of 
the related iteration processes) and not with intrinsic features of the Coulomb collision integral. 

This example of a successful application of our methodology to plasma oscillations supports 
the hope of further unraveling these and more entangled problems of plasma-echo phenomena. 
In particular, in the case of boundary conditions with two superimposed frequencies 



E(Q, t) = Eq cos(to>it) cos(a;2^) 
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3 For instance, in low-collision plasma. 



two non-damping waves with frequencies (u>\ — 0J2) and (uj\ + 0J2) arise in collisionless plasma. 
That is a very simple imitation of plasma echo effect (cf. description in ||) for both transversal 
and longitudinal oscillations. 
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